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Abstract. The Quantum Information Theory is a reach source of fascinating 
problems in Linear and Multihnear Algebra. In this paper we discuss one of 
such problems, namely the Distillation Problem. 
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^vj ' Let p^ , fc = 1, 2, . . . , m, be the critical Werner state in a bipartite dj. X d, 

^sj ' quantum system, i.e., the one that separates the 1-distillable Werner states 

from those that are 1-indistillable. We propose a new conjecture (GDC) as- 
serting that the tensor product of p^ is 1-indistillable. This is much stronger 
,,15 ' than the familiar conjecture saying that a single critical Werner state is in- 

^1^' distillable. We prove that GDC is true for arbitrary m provided that d^. > 2 

for at most one index k. We reformulate GDC as an intriguing inequality 
f*^ 1 for four arbitrary complex hypermatrices of type di X ■ ■ ■ X dm ■ This hyper- 

J^ ' matrix inequality is just the special case n = 2 of a more general conjecture 

'~{ ' (CBS conjecture) for 2n arbitrary complex hypermatrices of the same type. 

^^^ Surprisingly, the case n = 1 turns out to be quite interesting as it provides hy- 

permatrix generalization of the classical Lagrange identity. We also formulate 
the integral version of the CBS conjecture and derive the integral version of 
Ch ' the hypermatrix Lagrange identity. 
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10 ' 1- Introduction 

O, 

C^ \ There are two kinds of bipartite entangled quantum states. First, the states 

which satisfy the Peres-Horodecki criterion of separability (i.e., having positive 

semidefinite partial transpose) are known as PPT entangled states. Second, those 

which violate this criterion are known as NPT entangled states. The indistillable 

/\ ' entangled states are also known as bound entangled. (See the next section for 

>^ . the definition of distillability and 1-distillability.) The problem of existence of 

bound NPT entangled states was raised more than ten years ago [64 and is still 

open. There are several papers [3l [5l 13 [9l [13l |4] where evidence (numerical or 

theoretical) is provided for the existence of such states. Several researchers have 

proposed the conjecture 131 IS] that such states exist. We refer to it as the 

Distillability Conjecture (DC). This conjecture has been reduced [5] to the case of 

Werner states. More precisely, for a fixed bipartite d x d quantum system, it is 

known that bound NPT entangled states exist if and only if (iff) a Werner state 

with the same properties exists. 
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2 DRAGOMIR Z. BOKOVIC 

The non-normalized Werner states in a dxd quantum system can be parametrized 
as p^ (t) = I —tF, —1 < t < 1, where F is the usual flip operator (see its definition 
in the next section). They are separable for t < 1/d and entangled for t > 1/d 
[mil]. They are 1-distillable for t > 1/2 and 1-indistillable for i < 1/2 [HEIII]. 
Thus the state p^{l/2) separates the 1-distillable Werner states from those that 
are 1-indistillable. For this reason we refer to it as the critical Werner state. It 
has been shown that if the critical Werner state is indistillable, then this is also the 
case for all 1-indistillable Werner states [3]. 

In Section [3] we propose much stronger conjecture which asserts that the tensor 
product of several critical Werner states p^ — /3^(l/2), k = 1, . . . ,m, acting in 
a dk X dk system (with arbitrary finite dimensions dk) is indistillable. We refer 
to it as the Generalized Distillability Conjecture (GDC). We recall that there is 
an example fTTl of a distillable state ppyr <8) ^^^(1/2) for two pairs of qutrits (i.e., 
in the case rfi = ^2 = 3), where ppyr is a particular PPT entangled state and, of 
course, p^(l/2) is the critical Werner state. In view of this example it seems rather 
foolhardy to propose a conjecture like GDC. However, subsequent reformulation of 
GDC and a further sweeping generalization, as well as analysis and resolution of 
some special cases led us to this proposal. 

In Section|4]we reformulate GDC as an inequality for four complex hypermatrices 
of type di X • • • x dm- Let A^d denote the space of such hypermatrices, where 
d = {di,...,dm)- We embed this GDC inequality into an infinite collection of 
similar inequalities 

(1.1) $^"'(a;(i\...,x("),u(i),...,u("))>0, m,7i>l, 

where x^*^' , u^*"'' e Aid are arbitrary. The function $ ^ is defined by the formula 

$|j")(x(i),...,x("),u(i),.. .,«(")) 

where /„ = {l,2,...,m} and i = (ii, 12, . . . ,im) and j = (ji, J2, ■ • ■ , jm) with ik 
and jk running through {1, 2, . . . , d^}. 

We refer to the inequalities Eq. (|l.ll) as the generalized Cauchy-Bunyakovsky- 
Schwarz (CBS) inequalites, and we conjecture that all of them are valid (CBS 
conjecture). The case n = 2 corresponds to GDC. In the very first case, m = n = 1, 
the inequality Eq. (|l.ip is just the classical CBS inequality. 

In Section El we discuss some basic properties of the function $^" : M^ — > R. 
In particular, we introduce natural actions of the unitary groups \J{dk) and the 
general linear group GL„(C) on A^^" and we show that $j ' remains invariant 
under these actions. Another property shows that the cases where some dfc = 1 can 
be eliminated, i.e., the CBS conjecture can be reduced to the case where all d^ > 1. 

In Section [S] we prove that the inequalities (II. ip are true when m = I. We also 
show that they are true when n = 2 and at most one dk > 2. 

In Section [7| we prove that the inequalities (|l.ip are true when n = 1. Moreover, 
in that case we are able to express the function <I>^ {x,u) as a sum of squares of 
real valued polynomials. This gives the hypermatrix generalization of the classical 
Lagrange identity, see Eq. (|7.6p , valid over any commutative ring A. 
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As a concrete example, we write this identity in the case m = 2, i.e., when x and 
u are ordinary di x d2 matrices with entries in A: 

^i=l j=l J \k=l 1=1 / j=l 1=1 \i=l / 

2 
^ _ \ / d-, d^ \ 

i=i k=i \i=i 

= E 

1 < i < A: < di 
l<j<l<d2 

In Section [5] we formulate the integral version of the CBS conjecture and derive 
the integral version of the hypermatrix Lagrange identity. 

Finally, in Section [9] we summarize our results and conjectures. 

The symbols *, T and f denote the complex conjugation, the transposition and 
the adjoint, respectively. For any positive integer d we set I^ = {1, 2, . . . , d}. 

2. Preliminaries 

We consider a quantum system consisting of two parties, A and B (Alice and 
Bob), sharing a pair of particles. We denote by H = V.^ (g) H^ the Hilbert space 
for this pair. We assume that y,"^ and y,^ have the same finite dimension, d > 2. 
A product state is a tensor product n = p^ (^ p^ oi the states p^ and p^ of the first 
and second party, respectively. 

We remind the reader that a state or density operator, say p^, is a positive 
semidefinite operator T-l^ -> V.^ with unit trace, tr p"^ = 1. We often work with 
non-normalized states which are just nonzero positive semidefinite operators. A 
pure state \ijj) G ?^ is a unit vector; two such states are considered the same if they 
differ only by a phase factor. The density operator of the pure state \ip) is the 
1-dimensional projector p = \'4'){tp\; it is independent of the choice of the phase 
factor of \ip). We also refer to this p as a pure state. 

Any pure state \tp) € T-L can be written uniquely as \4') = ^Cy|i,j), where 
[cij] is a square matrix and we write \i,j) for |z) (g) |j). The Schmidt rank of |'0) is 
defined as the rank of the matrix [cij]. While this matrix depends on the choice 
of orthonormal (o.n.) bases of V.^ and V.^ , its rank is independent of this choice. 
One can choose these o.n. bases so that Cij = for i y^ j and cu > for all i. Then 
the Schmidt rank of \tp) is just the number of indexes i such that Ca > 0. The pure 
states lip) of Schmidt rank 2 play an important role in this paper (see the definition 
of 1-distillability below). 

We denote by p^ the partial transpose (T ^ l)(p) of a density operator p onH, 
where the transposition map T is computed with respect to some o.n. basis of H^. 

We shall use the following basic definitions. 

Definition 2.1. Let p be a density operator onV, — H^ ® Ti^ . 

(i) p is separable ij it can he written as a finite convex linear combination p = 
y^ AiTTi (Ai > 0, ^ Ai = 1) of product states tt^. It is entangled if it is not separable. 

(ii) p is PPT if its partial transpose p^ is positive semidefinite. It is NPT if it 
is not PPT. 
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(Hi) p is 1-distillable if {ip\p^\ip) < for some \tp) ^ H of Schmidt rank 2. It is 
1 -indistillable if it is not 1-distillable. 

(iv) p is distillable if p®"^ is 1-distillable for some m > 1. It is indistillable if it 
is not distillable. 

For more background information on entanglement, separability and distillability 
of bipartite quantum states we refer the reader to one of the papers [21 |3l El 13 [l2] 
or the recent book [1]. It is well-known that the separable states are PPT [TO] , 
It is also known that the distillable states are NPT ^. The question whether all 
NPT states are distillable has its origin in [7]. We shall refer to this question as 
the Distillability Problem. The answer to this problem is affirmative for d = 2, but 
widely believed to be negative for d > 3. Formally, the following conjecture, to 
which we refer as the Distillability Conjecture (DC), has been proposed in [3l Sec. 
II] (see also [i[S]). 

Conjecture 2.2. There exist indistillable bipartite NPT states. 

Let us fix an o.n. basis \i), i e Id of V.^, and an o.n. basis of V.^ for which we 
use the same notation. Due to the context, no confusion should arise. 
After fixing these bases, we can define the flip operator F : H ^ H hy 

where the indexes i,j run through Id, and we use the common abbreviations \x, y) — 
\x) (8) \y) and {x,y\ = {x\ ® {y\. The (non-normalized) Werner states on "H can be 
parametrized as 

(2.1) p^(f) = l-tF, -l<t<l. 

Let |(^) S H be the maximally entangled (pure) state given by 



W) = 



1 






\l,l) 



Its density matrix is the projector 

The partial transpose of p^ {t) is a^ {t) = 1 - tdP. 

The following facts about the Werner states are well-known. 

Proposition 2.3. The Werner states p^{t) are: 

(a) separable iff —1 < t < 1/d; 

(b) NPT iff 1/d <t< 1; 

(c) 1-distillable iff 1/2 <t<l. 

For (a) see (TSl [S], for (b) and (c) see [21 |H]. Since p^{l/2) separates the 
1-distillable Werner states from the 1-indistillable ones, we shall refer to it as the 
critical Werner state. The importance of Werner states for the distillability problem 
for bipartite states was first established in [6i. 

Proposition 2.4. For a fixed d x d quantum system, DC is true iff there exists an 
indistillable NPT Werner state. 
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We end this section with two examples. 

First, we verify that the Werner states p^{t) are 1-distillable for 1/2 < i < 1. 
We take \ip) = |1, 1) + |2,2), which has Schmidt rank 2. Then {iplip) = 2/Vd and 
PIV') = {2/Vd)\ip). Hence 

{i^W'^itm = (V'll-idPlV-) 

= 2 - 2tVd{Tp\ip) 
= 2(1 - 2t) < 0. 

In our second example we consider the Hermitian operators 1 — tdP where t is 
a real parameter. Since its eigenvalues are 1 and 1 — td, it is positive semi-definite 
iff i < 1/d. These states are known in the literature as isotropic states (see e.g. 
^^). The partial transpose of a separable state is also a separable state. Hence, the 
states a^{t) = 1-tdP = p^{t)^, -1 < i < 1/d, are separable. We shah prove that 
the states 1 — tdP, t < — 1, are 1-distillable. According to Definition 12. 1[ we have 
to show that {ip\l — tF\ip) < ioT some vector \ip) of Schmidt rank 2. The choice 
IV') = |1, 2) — |2, 1) works in all these cases. Indeed, we have {l — tF)\ip) = {l + t)\tp) 
and so 

(V'll - tF\xl:) = (1 + t)\\iPf = 2(1 + t) <0. 

3. Generalized distillability conjecture 



We now assume that Alice and Bob share m pairs of particles and denote by 



Hk = T-l-k ^ ^f the Hilbert space for the kth pair, k G Im = {1,2, .. . ,m}. We 



assume that the Hilbert spaces Ti^ and Tij^ have the same finite dimension, which 
we denote by dk ■ Although the case dk = 1 is not of interest for the distillability 
problem, we shall not exclude it for the sake of completness. The Hilbert space for 
the whole system of m pairs of particles is H — H^ (8) H^ , where 

We refer to d = (di, . . . ,dm) as the dimension vector of this composite quantum 
system. 

We fix an o.n. basis \ik), ik = l,---,rffe of H^, and an o.n. basis of "Hf 
for which we use the same notation. Let Fk : Hk ^ "Hk be the fiip operator 
and pY {t) = 1 ~ tFk the (non-normalized) Werner state. Let \(pk) € Hk be the 
maximally entangled (pure) state and Pk = \ipk){^k\ its density matrix. 

We assume that, for each k, the kth pair of particles shared by Alice and Bob is 
in some Werner state pY (tk) with tk < 1/2. We set t = ti, . . . , i^ and 

P^{t)=pY{tl)^pY{t2)^---^P^{t^). 

By using the same argument as in the proof of [3l Lemma 4] one can easily prove 
the following lemma. 

Lemma 3.1. If p^ {t) is 1-indistillahle and t' — t'^, . . . ,t'^ is such that t'^ < ti for 
each i, then p^(t') is also 1-indistillahle. 

In view of this lemma, we assume from now on that all pairs of particles shared 
by Alice and Bob are in the critical Werner states p^ — pY i^/'^)^ 1 < k < m, and 
we set 

p^^pY^pY^-'-^p"^- 
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We shall refer to p^ as a generalized critical Werner state. 

We are now ready to state our Generalized Distillability Conjecture (GDC). 



Conjecture 3.2. All generalized critical Werner states pY^ ^ ); 'ti > 1, 



Note that since m as well as the dk are arbitrary, if we replace "1-indistillable" 
with "indistillable" in GDC, we obtain an equivalent conjecture. In particular, 
GDC implies DC. 

We note that GDC is valid in the case m = 1 because the critical Werner states 
are 1-indistillable (see Prposition l2.3|) . A direct proof of a more general result will 
be given in Section [6l 

The partial transpose of p^ is 

(3.1) a^' = (1 - dlPl/2) ® (1 - d2P2/2) ® • • • ® (1 - draPml'2). 

Hence, GDC can be restated as saying that the inequality 

(3.2) (^|a^|V)>0 

is valid for all |V') G H of Schmidt rank at most two. 

There is a conjecture [21 [3j EJ [11] which asserts that the critical Werner states 
p^ {1/2) are indistillable. Equivalently, it asserts that the tensor power /9^(l/2)'*™ 
is 1-indistillable for all m>\. This conjecture is a special case of GDC, which can 
be seen by setting all dk equal to d. 

4. Generalization of GDC and the CBS inequality 

Our objective in this section is to reformulate GDC as an inequality for complex 
hypermatrices, and to embed this particular inequality in an infinite sequence of 
conjectural hypermatrix inequalities. 

Our first job is to find an explicit expression for the inequality p.2|) in terms of 
the components of jV'). We have |f/') — |^i) + |^2) where jV'i) and IV'2) are product 
vectors. Thus 

IV'l) = '^Xi^^,,,.i^Uj^^,,,^j^\ii,ji,...,i^,jm), 
ij 

IV'2) = ^2/n,...,i„Wji,...,i,„|«l,jl,---,im,jm), 

where i = ii, . . . , im] j — ji, ■ ■ ■ ,jm and, for each s, is and js run through the set 

Id,- 

We can also write \'4>i) — \x) (E) \u) and 1-02) — \y) <E} \v), where e.g. 

i 

Thus X = (xij_...^i,^) is an ?7i-dimensional complex matrix (or tensor with m in- 
dexes). For simplicity, we shall refer to such objects as hypermatrices. Let Md 
denote the space of complex hypermatrices of type di x • • • x dm,. Formally they 
can be identified with the space of complex valued maps on the Cartesian product 
-^d = ^di X • • • X Id,„. Finally we introduce the function $d : A^^ — > R by setting 

^d{x,y,u,v) = (V'lcr^lV')- 
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We have to find explicit expressions for each of the four terms ('0i|o'd^l'0j)i *i J G 
{1, 2}. We shall derive the one for i = 1 and j = 2. Let also p — pi, . . . ,pm and 
q — qi,...,qm where Ps and q^ will run through /^^ . We start with 

{iljl,..., irnJmW^\Pl,qi, ■ ■ ■ ^Pm^Qm) 

By expanding ()3.ip we obtain the formula 






E 



-1 



2iH hz„ 



0-21,1 (X) 0-22,2 <i$ ■■■<>$ rjz 



where z = (zi, . . . , z„j) and cro.^ = 1, (Ji,k = d^Pk- Using this we find that 



{^l\<\^2) 



QC/„ 



\Q\ 



ER i:'< 



, im'^jl,---,jmyPl '■■■■?!" ^91 '■■■.?™' 



where jQI denotes the cardinality of Q and the summation in S' is subject to the 
constraints 

Ps = is, qs = js for s G /,„ \ Q; j^ = is, Qs = Ps for s G Q. 

For instance, if to = 2 and Q = {2} the constraints give pi — ii, qi = ji, J2 — 12, 
Q2 =^ P2 and the sum E' can be written as 



iljl \ J2 

In general, E' can be written as 



\ P2 / 






^Jl,---Jn 



/ , yii,...,im^ji,...,j„ 

Kis,js=iB; s£Q 



While performing the last two summations one has to set first js ~ is for each 
s € Q and then to sum over all is, s € Q. Explicitly, using the Kronecker deltas, 
we have 

is,js=is;s£Q isJelsSiQ \q&Q j 

Similar three formulae are valid for (^i|o-^|'0j), (*,j) 7^ (1;2). By using all four of 
them, we obtain that 



(4.1) <I>d(x,y,M,i;)= E ( — ) E 



-1\ 



\Q\ 






Let us give two concrete examples. If to = 1 and d = di, then 



1 



i,i=l 



^(a;iUi + yiWi) 
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Writing x, y, u, v as row vectors and using the Frobenius matrix norm, this can be 
expressed as 

'^{d){x,y,u,v) = \\x'^u + y'^v\\'^ - - \xu'^ + yv'^\ . 

If TO = 2 then x,y,u,v are ordinary di x d2 complex matrices and 



'^idi,d2){x,y,u,v) = ^ \x^JUkl + yijVkil'^ - -Y^ 



i,j,k,l 



j,i 



Y^ixijUii +yijVii) 






J 



{xijUkj +yijVkj) 



/ X ^ijUij + yijVij) 



«j 



where the indexes i, k run through Id^ and j, / through /^j . Using the tensor product 
of matrices, this can be written in a more compact form 

*(di,d2)(a;,y,'«,w) = \\x®u + y®vf - -\\x'^u + y^vf 

--|lua;^ + wy^||^ + -|tr [x^ u + y'^ v)\^ . 

In the special case di = d2 = d this formula was derived in j4j . 
From Eq. (|3.2p we see that GDC can be restated as follows. 

Conjecture 4.1. The inequality ^d{x, y, u,v) > is valid for all dimension vectors 
d = (di, . . . , dm), TO > 1, and all complex hypermatrices x, y,u,v (z Add- 

We shall now embed this inequality in an infinite sequence of hypermatrix in- 
equalities. The sum XiUj + j/iUj occupies a prominent part in the definition of the 
function $d : Ai^ ^- R given by Eq. (|4.ip . Clearly, it is rather unnatural for 
this sum to have only two terms. Therefore we introduce a more general function 



$ 



(») 



:M2" 



R such that <i)d = <i> 



— ff>^ . It is defined by a similar formula where the 



sum XiUj + j/ifj is replaced by one which has n summands, and at the same time the 
fraction —1/2 is replaced with — 1/n. As we will see in Section[5l this new function 
shares with <I>d several properties observed in [4J in a special case, including the 
conjectural property of nonnegativity. The new function is defined as follows: 



(4.2) 



M(M) 



^'d'i^ 



,x("),u(i) 



.,«(")) 



Ef^ 



\Q\ 



E E-i 



(k) (fc) 

u- 



iqJq=iq;qeQ 



k=l 



E 

where d = (di , . . . , dm) and x^*^' , m*^*^^ & M.d are arbitrary hypermatrices 
We now state the Cauchy-Bunyakovsky-Schwarz (CBS) conjecture. 

Conjecture 4.2. The inequality 



is valid for all dimension vectors d 
hypermatrices x^^\u^^' G A^d- 



,,a:("),«(l). 



(di,. 



...,w("))>0 

dm), all m,n > 1, and all complex 



We will see later that this conjecture is true (and nontrivial) in the two boundary 
cases TO = 1 and n = 1, see Sections[n]and[7]respectively. In the case to, = n = 1 this 
is indeed the classical Cauchy-Bunyakovsky-Schwarz inequality, which explains our 
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name for this conjecture. We also point out that the inequahty in Conjecture 14.11 
is just the case n = 2 of the CBS conjecture. Thus the CBS conjecture generahzes 
both the GDC and the CBS inequahty. 



5. Some properties of the function ^^j"-^ 

We generahze here the properties of ^d noted in [4 in a special case. 

5.1. Generalization of matrix transposition. The transposition operation on 
ordinary matrices admits a wide generalization to hypermatrices. Let us fix a 
permutation tt of the set /„ = {1,2,..., m}. For a given dimension vector d = 
(di, . . . , dm) we set d' = {d'l, . . . , d'^), where d'f, — d^-ife, k £ Im- For x G Ma we 
define x' £ M.^' by 

This is justified since i'^ G Id' = Id _i implies that i^j, G Id^- Using the abbrevi- 
ation a; = (^x^''^) , we have 

,x'("\u'W,. ..,«'(")) 
= <i>^"^(x(i),...,x("),Ji),...,w(")) 
for all x^^^u'^^^ G Md- 



(5.1) 



<"^(-'^^ 



5.2. Elimination of ds = 1. We show that the ds which are equal to 1 can be elim- 
inated. Assume that ds ~ 1 for some s G /,„. Set d' = (rfi, . . . , ds-i,ds+i, . . . , dm)- 
For each x G M.d we define x' G A^d' by setting 



.,Ia_l,Ja + l, 



f'ii,...,is_i,l,is^i,...,i„ 



Note that this x' is different from the one used in the previous subsection. 

Proposition 5.1. Using the above assumption and the abbreviation x = {x^^') , 
the equality 

<i>(")(x«, . . . ,x("),««, . . . ,u(")) = !i^$('^)(x'^^\ . . . ,x'^"\«'^^\ . . . ,«'("^) 

n 

is valid for any x^^\ u^^' G A^d- 

Proof. Note that 



(5.2) 

where 
(5.3) 



/ 1 \ l'3l 



$(;')(x(i),...,a;("), «(!),...,«(")) 



E 



ip,jp\P£lm\Q 



E E-l 



,Wy(fc) 



We set Im — Im \ {«}. Let Q C /„ be a subset of cardinality \Q\ = k + 1 such that 

Q' 



s G Q, and set Q' ~ Q\ {s}. Because d^ = 1 we have $q ' = ^q' ■ Hence 



-1 



fc+i 



$1."' 



-i^%w^!i^.(^ziVcogV 
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It remains to observe that 



(")/^'(l) „'(») „,'(!) 



ci>^:^(.' 



X ,u 



,u'(")) 



'C/' 



IQ'I 



E (-)■"■<(-«,. ..,x("),uW,. ..,«(")). 



We single out the foUowing special case. 

Corollary 5.2. If d — (di, . . . ,dm), m > 1, and some dg = 1, then for any 
x,u £ A4d we have $^ {x,u) — 0. 

5.3. U((is)-invariance of $j . We introduce an action of the unitary group \J(ds), 
1 < s < m, on the space A^d- For U S \J{ds) and x G A^d we define [/'"■'a; G A^d 
by the formula 

To simplify the notation, we set U*^^'^x — {U*y'x. 
Proposition 5.3. $^' is \J{ds) -invariant, i.e., 

/or aZ/ a;W,u('=) G Md and all U G U(4). 

Proof. In view of Eq. (|5.2p , it suffices to show that for each Q C /,„ we have 

$[j") ([/(^)x(i) , . . . , C/(^)x(") , f/*(^)u(i) , . . . , f/*(^)u(")) 

= <i>J'(xW,...,x("\w(i),...,u(")), 

where $q is defined as in Eq. (|5.3p . 

If s G Q this is easy to verify. Indeed we have 



2 


^ (ty-'a:ji(ty--ujj = 2. Z^^»=. 


»',f^»,.j- 2;«i 


■■■,Js.---.*m%'l.---j'a.---J" 


is -.3 s— is "^'s^^'s ^ '^s 


/ 




■,ii,---,Jm%'l,---,i^,---,im 


is.js—'i-s 




li s ^ Q then we consider the sum 


2 


1. 


^ ^([/(^)xW)i(C/*(^)tiW)j 






»s J's 


iq,jq=iq;q&Q k=l 



\uzum 



where Z is the matrix with entries 



E E -£.... 



,(fc) 



■,*Tn Jl,--- Js,--- Jj 



,; ■i'sJ's&Id, 



iq,Jq=iq;<l&Q k=l 
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Since U is unitary we have 

2 



\uzu^f = \\zf = Y, 



which completes the proof. H 

The actions of the unitary groups U(ds) for s G /m on Aid pairwise commute, 
and so we obtain an action of the direct product U(d) = U((ii) x • • • x U(dm) on 
A^d- Clearly, ^J^ is also U(d)-invariant. 

5.4. GL„-invariance of $|j"\ Let A = [ap.q] £ GL„(C) and (A^)-i = [(3p.q]. 
Proposition 5.4. Using this notation, we have 



\fe=l k=l fe=l fe=l 

(")('i^(l) -r-(") -7,(1) ,/(") 



for all x^^\u^^^ e Md- 

The proof given in [U Section 4] for a special case extends easily to the general 
case. 

6. Some special cases of the CBS conjecture 

We start with a simple direct proof of the CBS conjecture in the case m = 1. In 
particular, this gives an independent proof of GDC when to = 1. 

Proposition 6.1. The CBS conjecture is true when m — 1. 

Proof. We set di — d and represent the hypermatrices a;*^''\ u*^*^-* S A^(d) as row 
vectors. The CBS conjecture takes the simple form: |jX|p > (l/7i)|trXp, where 

is a dx d matrix of rank r < n. By Schur's triangularization theorem, we can choose 
B e U(d) such that Y — BXB^ is upper triangular. Then the diagonal entries 
Ai, . . . , Arf oiY are its eigenvalues and, by the same theorem, wc may assume that 
Aj: = for r < A; < d. Since \\X\\ = ||y|| and tr X = tr y, our inequality becomes 
\\yf > il/n)\Xi + ■■■ + Xr\^. Since \\Yf > \Xi\^ + ■ ■ ■ + \Xr\^ and r < n, the last 
inequality is a consequence of the well known inequality 

Ml H 1- Mr ^ /^Mi H ^A*r^ 



where mij • ■ • j Mr are nonnegative real numbers. ■ 

In the case n = 2 we can extend Proposition IG . 1 1 by making to arbitrary provided 
that di > 2 for at most one index i. However in this case we rely on the equivalence 
of GDC and the case n = 2 of the CBS conjecture. 

Proposition 6.2. Let pi, i — 1,2, be a bipartite state acting on Hi = 'Hf®'Hf . If 
Ml is separable and pi 1-indistillable, then pi ® p^ is 1-indistillable. Consequently, 
CDC is valid for arbitrary m provided that di > 2 for at most one index i. 
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Proof. Let Ci = (T (g) l)pi, i — 1,2, be the partial transpose of pi. Since pi 
is separable, cti is positive semidefinite and separable. Let {ip) S Hi (8) ^2 be any 
vector with Schmidt rank < 2 (with respect to the partition A\B). We have to show 
that {ip\<7i (E) (T2\tp) > 0. If tr 1 is the first partial trace function (i.e., corresponding 
to Hi), then 

ii^Wi (E) a2\ij) =tr (|i/^)(V'|cri«)fT2)=tr [tr i (|7/')(V|cri) (72)] • 

We now use an argument from [5l Section III B]. Since ui is separable, we can write 
it as CTi = X]iCi|aj,6j)(ai,6j|, Ci > 0. Then 

i 

where |^i) = (a^, 6i|'0) is a state acting on 'H2- Next we have 

{il>\ai (g)cr2|^> = ^Ci{'il:i\(j2\il)t). 

i 

As the local projections cannot increase the Schmidt rank, each \'ipi) has Schmidt 
rank < 2. Since p2 is 1-indistillable, each summand in the above sum is nonnegative. 
This concludes the proof of the first assertion. 

The second follows from the first and the Eq. (|5.ip . ■ 

The proof of GDC can be reduced to the case where 3 < di < ^2 < • • • < dm- 
This is the objective of the next proposition. 

Proposition 6.3. // Conjecture \4-l\ is true for all d = (di, . . . , dm), m > 2, 
satisfying 3 < di < d2 < ■ ■ ■ < d„i, then it is true in general (i.e., without any 
restrictions on the dimension vector d). 

Proof. Let d = (di, . . . , dm) be arbitrary and let p^ be the corresponding gen- 
eralized critical Werner state. We have to prove that ^d{x,y,u,v) > for all 
a;, y,u,v G ^Ad■ Because of Eq. (|5.ip we may assume that di < d2 < • • • < d„j. By 
Proposition 15. li we can further assume that 2 < di . If 3 < di then the assertion 
follows from the hypothesis. 

Otherwise we have 2 — di = ■ ■ ■ ^ dk for some k < m. The bipartite state p' = 
pY E) ■ ■ -^ pY is separable, lik — m then p^ = p', and so it is separable and hence 



IS 



1-indistillable. Finally, let k < m. Then the bipartite state p" = Pk+i ® • • • ® p, 
1-indistillable by hypothesis or (if m- = A: -I- 1) by Proposition 12.11 P^^rt (c). Hence, 
p^ — p' E) p" is 1-indistillable by Proposition 16.21 ■ 



7. Hypermatrix Lagrange identities 

It will be shown in this section that the CBS conjecture is true in the case n 
Indeed, this is an immediate consequence of Theorem 17.21 below. 
From Eq. (|4.2|) we see that 






^T 



xm 






where d = (di,...,dm) and a;,M G A^d- We shall prove that the inequalities 
<i>|j^^(x,u) > indeed hold. 
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We write in detail the first two cases. For m = 1 we set d = di and obtain 

2 



<fS(-,-)-Ei-^i'-E 



U-j 



i=l 



j=i 



E 
fc=i 



XkUk 



The hypermatrices x and u are just vectors in C . 
For TO = 2 we have 



/ , •^n,J2^n j2 



E 



/ , •'^ilM'^jl, 



+ 



/ ^ Xi^^i^Ui-^^' 



where ii and j'l run through Id^ and «2 and j2 through 7^2 • In this case the 
hypermatrices x and u are ordinary di x d2 complex matrices. 

The inequality ^,l{x,u) > is the well known Cauchy-Bunyakovsky-Schwarz 
inequality, which follows trivially from the classical Lagrange identity 

2 



EN^-Ei 



i=l 



J = l 



E 
fc=i 



XkUk 



E 

2j'; i<3 



,*|2 



.(1) 



The inequality $L ^ \(a;, m) > follows trivially from the following matrix ana- 
log of Lagrange identity: 



.(1 



$7 . Ax,u) 

{di,d2)^ ' ' 



/ y \-^il,i2^ji. 



J2 ^iuJ2^Ji,i2 



* 1^ * 1 2 

^jl,i2'"ii,j2 ^^ ^ii J2"ii,i2 I ' 



ij;i<j 



where we recall that i ~ ii,i2 and j ~ ii,.72j and the inequality i < j means that 
ii < ji and Z2 < j2. 

There is also an analog of Lagrange identity for any m. To state and prove 
these identities in general we need some additional notation. Let i ~ ii, . . . ,im and 
j — ji, . . . , jm be two sequences of indexes with is^js G Id^ for each s. For arbitrary 
Q (- Ijn we define two modified sequences of indexes i*^'J and j'^^' of length to. as 
follows: 



iQJ = / Jp ifpeQ; 



ip otherwise; 



JP 



iipeQ; 
otherwise. 



(7.1) 

We also define 

(7.2) ai-j{x,u) = y^ {-ly^^XjQ.jUp.i. 

For TO, = 1, 2 we have: 

^il,i2;jl,J2\-^T^ ) '^*l'*2^Ji J2 Xij^J^U^ 

Lemma 7.1. Let i = ii, . . . , i^ and j — ji, . . . ,jm be two sequences of indexes with 
isijs € Id, for each s. For any Q Q Im and any x,u E Md we have 

a\'-j\x,u) = (-l)l'5lcri;j(a;,M), 



"^1,^2 •"■'1 J2"-J"i,i2 ^jl,i2'^il,J2 ^ ^h,J2''^ix.i2- 
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where i' = i'^'J and j' = j*^''. 

Proof. We first observe that ii Q, RC /,„ and 

^ ^ : J J ; 

then 

•// _ iQARj •II _ •QAR,\ 

^ ^ 5 J J ; 

where QAR is the symmetric difference of Q and R. Hence 

ai'.y{x,u) = ^ (--l)'^'xi,Hj'Mj,H.i' 
_R,C/,„ 

— 2_^ ("1) a^iQAH.jUjQAH,, 

= (-l)l'3lai;j(a;,u), 

where we used the congruence |QA_R| = |Q| + \R\ (mod 2) and the fact that when 
R runs through all subsets of Im so does QAR. ■ 

In particular, note that if ig = js for some s then cri;j(a;, u) = 0. 

We can now state our hypermatrix Lagrange identity in full generality. 

Theorem 7.2. For any dimension vector d = {di, . . . ,dm), rn > \, and any 
hypermatrices x,u d Ai^ we have 

(7.3) <i>W(^,^) = _L^|^..(^,„*)|2= J2 ki;j(x,u*)^ 

where i < j means that is < js for all s. 

Proof. To prove the second equality, let 51 be the set of all pairs (i; j) such that 
is y^ js for all s, and let 51^ be its subset consisting of all pairs (i;j) such that 
is < js for all s. Since ai-j{x,u*) = if is = js for some s and, by Lemma 1 7. 11 

|CTiQ,j.jQ,i(x, U*)| = |cri;j(x, u*)| 

for all Q Q Im, we have 

i.j (i;j)en 

= 51 X! |CTiQ,J;jQ.i(x,U*)| 

(i;j)en# QC7„, 

= 2™ Y. I^id(^' «*)!'' 
(i;j)6n# 

= 2™ Y. I^i;j(^'^*)l'- 
The proof of the first equality is more involved. From Eq. (|7.2p we obtain that 
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Hence 

(7.4) 5]|aij(x,7.*)|'= Yl (-l)"^'''"E^i«'^"J«''<«.J"I«.^ 

From Eq. (|7.ip we see that, for a fixed s e /„, we have 

(«s,«s,js,js) ifsGi?\Q; 

{js,js,is,is) iiseQ\R; 

{js,is,js,is) liseQnR; 

{is,js,is,js) otherwise. 

By using this, it is not hard to deduce that 






(7.5) ^a;iQ,jUji?..iXj*fi,jUj*Q.i = ^ 



y xiQ.jUjH, 

is,js=is;s£QAR 



The correct interpretation of this equality requires that we omit the summation 
sign on the right-hand side if the set of indexes over which we sum is empty. For 
instance let m = 2 and Q = R= {I}. Then i'^J = i^^J = ji, ij and j'^'' = j-^^^-' = 
n, J2, and we have 

a^iQjMjH.ia^iRjUjQ.i — Z^ ^Jl ''2^11 ,J2^Ji ,12^*1 J2 

iJ 

This agrees with Eq. (|7.5p provided that 

/ XiQjUjfl.i = XiQ.jUji?..i (= Xj^.i^Ui^J^) ■ 

iB,JB=is;s£QAR 

According to our interpretation, this equality is indeed valid since QAR — 0. 

Since for a fixed S C I^ there are exactly 2™ choices for subsets Q,R '^ Im such 
that QAR = S,we deduce from Eqs. (fTI)) and (fTS)) that 

2 






5:Kj(x,.*)r=2™^ (-1)1^1 5: 

iJ SC/„, isJs;se/m\S 

This completes the proof of the theorem. ■ 

Observe that if some d^ = 1 then the inequality i < j cannot be satisfied and so 
Eq. (fTS)) gives ^^J^\x,u) = 0. This agrees with CoroUarv [5?2] 

The following purely algebraic version of the generalized Lagrange identity is 
valid in a more general context. 

Corollary 7.3. For any dimension vector d — (di, . . . ,dm), rn > 1, and any 
hypermatrices x, u of type di x ■ ■ ■ x dm with values in an arbitrary commutative 
ring, we have 

(7.6) ^(-1)1^1 J2 E ^i«J"j«' - E ^i;j(^'")'' 

QQIm. ipdp;peI,„\Q \«<!;9eQ / ij:i<j 

where we make use of the definitions in Eqs. {Ol) and (TW). 
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Proof. It is immediate from the theorem that the identity Eq. (|7.6p is vahd when 
the hypermatrices x,u are real- valued. Since this is a polynomial identity over the 
ring of integers, it is also valid when x and u are hypermatrices taking values in an 
arbitrary commutative ring. ■ 

8. Integral version of the CBS conjecture 



Let _L^(R™) be the space of complex- valued square-integrable functions ^ : 



R" 



C. We define the function $^^ : L2(R™)2" ^ R by the formula 



<^l:H^^'\...,^^-\v^'\...,v^-^) 



E 

0C/,„ 



n I 



\Q\ 



Sp,tp;p£lm\Q 



Sq.tg — Sq-^q^Q ^._-|^ 



i^'^KM^Ht) 



where ^C^), r/C"') e L^(R'") and s = si, . . . , s^ and t = ti, . . . , i,„. While computing 
the inner integrals (those within the absolute value signs) one has first to set tq = Sq 
for each q (z Q and then to integrate over each variable Sq, q G Q. The outer 
integrations are performed over each of the variables Sp and tp with p e Im\Q- 

In the case when the functions ^'*''' , i]^'^' are continuous with compact sup- 
port, one can approximate them by hypermatrices x^*'\u^*'' e A^d where d = 
{di, . . . ,dm)- After suitable scaling and by taking the limit as all dg — > oo, we 
obtain a version of the CBS conjecture in which all the sums are replaced by inte- 
grals. Since the continuous functions with compact support are dense in L^(R'"), 
the CBS conjecture implies that 

(8.1) <i>(;:)(eW,...,e("\r;W,...,r;("))>0, 

for arbitrary ^('=),7y('=) e L2(R'"). 

For instance li m — n — 2 then the inequality reads as follows: 



dsi ds2 dti dt2 

2 



(£.^^\suS2W^\siM)+S,^^\si,S2W\siM))dsi 

[^^'\si,S2WHh,S2) + i^^Hs,, S2W^\h, S2)) dS2 
(S,''^\si,S2W^\si,S2)+i^''\si,S2W^Hs^,S2))dsidS2 



ds2 dt2 

[ 

dsi dti 

2 

>0, 



where each integral is taken over the whole real line. The four expressions on the 
left-hand side are contributions of the four subsets Q = 0,{1},{2},{1,2} of/2, 
respectively. However, we want to point out that this particular inequality is only 
conjectural. 

Since we know that the CBS conjecture is true if m = 1 or n = 1, the inequality 
T]) is valid in these two cases. Explicitly, these inequalities are: 

$(")(^(l) f:{n) „(1) 



,e^"^v 



k=l 



e^'isWit) 



(k), 



,vn 



ds dt 



fc=i 



e^'=^(s)ry('=)(s)ds 



>0 
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and 



QC/„ '' ''Sp,tp 



P<^I^\Q 



?(s)^(t) 



Sg,tq=Sq; gGC 



>0. 



The latter inequality is also a consequence of the integral version of Eq. (I7.3P 
which we are going to write down explicitly. Recall that s — si, . . . , s™ and t ~ 
ii,...,im are sequences of real variables. For Q Q Im we define two modified 
sequences, s*^'* and t'^'^, of the same length m: 



tp if p G Q; 
Sp otherwise; 



Sp if p e Q; 
tp otherwise. 



We also define 



as;t{i,v)^ E (-i)"^'^!^^'*)^!^^'^ 



3C/„ 



Now the integral version of the hypermatrix Lagrange identity (17.31) can be stated 
as follows: 



^ "' ■JSr,,tr,;pGlm 



If 



ks;t(^,r?*)r 



Sp<tp;'pel,r 



where ^, r/ S L^(R'") are arbitrary. 

We can derive additional inequalities from the CBS conjecture by letting only 
some of the dimensions ds — >■ (X) and keeping the other fixed. As an example 
we consider the case where m = n = 2 and d = (^1,^2), where di ^ 2 while 
we let ^2 -> 00. By Proposition 16. 2[ we have ^d{x'^^\x''^\u''^\u''^^) > for aU 
j;(i)^2;(2)^ y(i)^^(2) g M.d- In this case we obtain the inequality 



i 

iE/(^r^(«)'?r^(^)+^f^(^)^r^(^))^^ 



ds dt 



>0, 



valid for arbitrary functions Q ,1]) ^ -^^(R-)j k — 1,2. (As di = 2, the indexes i 
and j run through /2 = {1,2}.) 

We give two concrete examples. 

First, let ^l''\t) = i(°.'=-i)/2 and 77^ ^(t) = t^'',k-i)/2 .^yi^gj.g < i < 1 and aU a,fc 
and bik are positive. We assume that all of these functions vanish for i < and for 
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i > 1. Then the above inequahty becomes 



4 V i 2y 



i,j,k,i ("*'= + ^»')(^Jfc + ^J') ~^^i ('^'fc + '^jOibik + bji) 



~' \f^ ^ik + bjk J I Tfe' '^*'= + ^'■'^ 



>0, 



where i, j, k, I run independently through {1, 2}. 

Second, let Q (t) — e^"-^*'*' and rji ' [t) — e~''*'=* where all aik and bik are 
positive. Then we obtain the inequality 

4V ^ ' -2V 



^(flife + aa){bjk + fcjO J , J. ; ^(fljfc + aji){bik + &j;) 

2 

>0. 






\/aik + b]k. J \ J j^^ V(Hk+bik 



i.j \ k V ^'^ J'" / \ i,k 



9. Conclusion 



The Distillation Conjecture (DC) asserts that there exist indistillable NPT (nec- 
essarily entangled) states in a bipartite dx d quantum system with d > 3. Although 
it is more than ten years old, it remains open for all d > 3. DC has been reduced to 
the case of non-normalized Werner states p^ (t) = 1 — tF, — 1 < i < 1, where F is 
the flip operator. For the affirmative answer to DC, it suffices to show (and it has 
been conjectured) that the critical Werner states p^ = p^(l/2) are indistillable. 
These critical Werner states are known to be f-indistillable and lie on the boundary 
between the f-distillable and the 1-indistillable Werner states. 

We propose much stronger Generalized Distillation Conjecture (GDC) which 
asserts that the tensor product of arbitrary critical Werner states p^ of type dkxdk, 
k = 1, . . . ,m, is 1-indistillable. It is easy to see that "l-indistillable" can be replaced 
here by "indistillable" . We recall from [TT that there is a distillable bipartite state 
ppyr (8) ^^^(1/2), where ppyr is a particular PPT entangled state of two qutrits 
and p^(l/2) the critical Werner state of two qutrits. In view of this example, 
GDC appears to be somewhat counter-intuitive as it asserts in particular that 
^^^(1/2) (E) p'^ (1/2) is indistillable. However, there is a strong numerical evidence 
that the latter state is indistillable. 

We reformulate the GDC as an inequality ^d{x,y,u,v) > 0, see Conjecture 
14. H where d — {di,...,dm) and x,y,u,v are arbitrary complex hypermatrices 
of type di X • • • x d„i. We are then led to propose a more general inequality 
*!^^(2:^^^ • • • ,a;("),u(i), • • • ,w(")) > 0, see Conjecture S^l which includes GDC as 
the special case n = 2. It also includes the classical Cauchy-Bunyakovsky-Schwarz 
(CBS) inequality as the case m = n = 1. For that reason we refer to the latter 
conjecture as the CBS conjecture. 

We have shown that the CBS conjecture is true in the two boundary cases 
m — 1 and n = 1. The case n = 1 is quite interesting since we were able to express 
the function $j (a;,u) explicitly as the sum of squares of real- valued polynomi- 
als. These polynomial identities provide generalizations of the classical Lagrange 
identity, which is obtained for m = 1. 
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As in the well-known case of the classical CBS inequality, the discrete inequalities 
in the CBS conjecture have their continous counterparts in which the hypermatri- 
ces are replaced with complex- valued square-integrable functions on R™, and the 
summation operations are replaced with suitable intagrals. We formulate explicitly 
these integral inequalities and point out some special cases including two concrete 
examples. We also give the integral version of the generalized Lagrange identities. 

We hope that our conjectures and results will stimulate further research on the 
original Distillation Conjecture. 
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